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An  asymmetric  dielectric  barrier  discharge  model  is  presented  for  real  gas  eight  species  air  chemistry  using  a 
self-consistent  multi-body  system  of  plasma,  dielectric  and  neutral  gas  modeled  together  to  predict  the 
electrodynamic  momentum  imparted  to  the  working  gas.  The  equations  governing  the  motion  of  charged  and 
neutral  species  are  solved  with  Poisson  equation  on  finite  elements  using  a  Galerkin  weak  formulation.  Two 
specific  geometries  were  used.  First  the  unsteady  electrodynamic  force  development  mechanism  is  studied  over  a 
flat  plate  due  to  charge  and  neutral  species  production  from  adjacent  air-like  N2/02  mixtures  in  a  radio 
frequency  driven  barrier  discharge.  The  solutions  show  a  dominance  of  CT  and  N2+  ions.  The  time  average  of  the 
streamwise  force  shows  mostly  acceleration  above  the  actuator  but  a  small  decelerating  force  downstream  of  the 
powered  electrode  inducing  a  fluctuation  in  the  temporal  evolution  of  streamwise  velocity.  This  decelerating 
component  needs  further  investigation  as  this  may  inherently  limit  the  performance  of  a  dbd  plasma  actuator. 
The  computed  electron  temperature  based  on  the  electric  field  resembles  the  visible  glow  that  is  commonly 
reported  in  the  literature.  In  another  case,  neutral  gas  flow  is  considered  over  a  5  deg  conical  crosssection  for 
demonstration  of  active  control  using  pulsed  direct  current  and  dielectric  barrier  discharges.  The  results 
demonstrate  the  need  for  utilization  of  multiple  actuators  for  controlling  a  largely  separated  flow. 
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Subscripts: 

Number  density,  cm'3 
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Species 

Nitrogen  molecule 

B 

Boltzmann 

Nitrogen  atom 

E 

Electron,  finite  element 

Positive  nitrogen  ion 

I 

Ion 

Oxygen  molecule 

N 
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I. 

Introduction 

ELECTRODYNAMIC  modifications  of  air  flowing  around  flat  plates  and  air  craft  wings  have  been  documented  both 
experimentally  and  numerically.1"3  These  are  primarily  surface  barrier  discharges  using  asymmetric  configurations 
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of  electrode  sets.  The  symmetric  configuration  of  electrode  sets  can  produce  flow  oscillation  similar  to  synthetic  jets. 
For  plasma  actuators  with  spatially  displaced  electrodes  the  airflow  inducement  is  close  to  the  dielectric  surface  with  a 
time-averaged  flow  velocity  measured  between  1-3  m/s.  These  actuators  exhibit  several  benefits  in  active  flow  control 
applications,  including  absence  of  moving  parts,  rapid  on-off  deployment  and  attractive  self-limiting  characteristics. 
The  discharge  may  be  characterized  as  transient  with  a  time  scale  of  the  order  of  a  few  microseconds  while  that  for  the 
bulk  gas  flow  is  in  milliseconds.  This  allows  the  fluid  to  respond  relatively  instantaneously  to  the  electrical  inputs. 
Experimental  study  and  numerical  predictions  on  the  effect  of  different  discharge  parameters  such  as  applied  voltage 
waveform,  distance  between  electrodes,  dielectric  thickness  and  permittivity  are  documented  by  various  researchers.4"5 
These  results  show  the  volumetric  body  force  increases  as  a  function  of  frequency  and  applied  voltage.  The  parametric 
relationship  is  not  clear  and  is  an  area  of  continued  investigation.  Also  the  nature  of  the  plasma  gas  interaction  is  not 
well  understood.  Acoustic  measurements  of  the  dielectric  barrier  discharge  (DBD)  plasma  actuator6  suggest  that 
compressibility  effects  may  contribute  in  the  momentum  coupling  between  plasma  and  neutral  gas.  Specific  modality  of 
discharge  whether  it  is  glow,  Townsend  or  streamers  will  help  determine  how  the  body  force  is  transported  from  charge 
to  neutral  species.  The  effect  of  this  body  force  in  controlling  three  dimensional  laminar  and  transitional  flow  structures 
that  have  been  studied  utilizing  both  phenomenological  and  first-principles  models. 

Phenomenological  and  circuit  based  models  have  shown  some  success  in  mimicking  overall  electrical 
characteristic  of  the  discharge  but  are  fundamentally  limited  from  inherent  static  and/or  predetermined  charge 
imposition.7  The  first  principle  models  suitable  for  capturing  detailed  dynamics  are  also  computationally  challenging. 
For  example,  the  surface  barrier  discharge  in  these  actuators  consists  of  repeated  burst  of  micro-discharges  which 
appears  as  an  optical  glow  due  to  nanosecond  timescales.  Resolution  of  a  typical  micro-discharge  driven  by  electron 
decay  requires  running  simulation  at  a  timescale  of  picoseconds  or  lower,  thus  imposes  practical  restriction  to  ever 
reaching  the  fluid  timescale  which  is  in  milliseconds.  Recent  development  of  a  loosely  coupled  fluid  plasma 
formulation5,8  has  shown  reasonable  success  in  describing  radio  frequency  (rf)  powered  electro-hydro-dynamic  (EHD) 
control  of  NACA  wing  stall,  laminar  separation  over  a  ramp,  boundary  layer  transition  on  a  flat  plate,  and  turbulent 
separation  over  a  wall-mounted  hump.  The  plasma  equations  are  based  on  a  self-consistent  simulation  of  charge  species 
continuity  and  electric  field  equations.  Once  the  reasonable  periodicity  is  achieved  in  plasma  characteristic,  the 
temporally  resolved  EHD  force  is  employed  into  the  fluid  equations  to  predict  the  effect  of  plasma  in  flow  actuation. 
Interestingly,  for  moderate  Reynolds  numbers  either  co-flow  or  counter-flow  pulsed  actuators  with  sufficiently  high 
frequency  (kHz)  were  found  to  mitigate  wing  stall.9  Unsteady  flow  actuation  with  a  duty  factor  seemed  to  perform  better 
than  continuous  operation  of  the  rf  plasma  actuator.  These  results  underscore  the  importance  of  accurately  resolving  the 
unsteady  body  force  and  understanding  its  transition  and  turbulence  enhancement  mechanisms  for  effective  plasma 
actuation. 

Numerical  simulations  of  an  asymmetrically  arranged  dbd  actuator  employing  multispecies  hydrodynamics5 
and  Monte-Carlo  (PIC-DSMC)  methods10  by  independent  research  groups  have  found  that  ionization  is  not  equal  during 
positive  and  negative  part  of  the  alternating  cycle  producing  a  net  force  in  one  direction.  Recently,  the  dc  corona 
discharges  and  ac  dielectric  barrier  discharges  have  been  investigated  under  a  range  of  operating  conditions.  Adjustment 
of  the  actuator  geometry,  dielectric  materials,  power  frequency  and  RMS  voltage  show  improvement  in  power 
transferred  to  the  neutral  gas  flow  by  ion-neutral  collisions.  Also  power  losses  due  to  inadequate  impedance  matching  of 
the  power  supply  to  the  actuator,  dielectric  heating,  and  power  required  to  maintain  the  atmospheric  pressure  plasma 
have  been  reduced  for  better  performance  of  actuator.4 

In  our  earlier  work,  we  have  shown  that  electron  deposition  downstream  of  the  overlap  region  of  the  electrodes 
results  in  formation  of  virtual  negative  electrode  which  always  attracts  the  ions.  The  force  on  the  charge  separation  is  in 
the  positive-x  direction  during  most  of  the  positive  half  of  the  rf  voltage  and  in  negative-x  direction  during  most  of 
negative  half  of  the  rf  voltage.  The  magnitude  of  the  former  force  is  higher  that  the  latter.  Thus,  domain  integration  of 
the  force  on  the  charge  separation  and  time  average  of  the  force  was  found  primarily  to  be  in  the  positive-x  and 
negative-y  direction.  We  solved  the  equations  governing  dynamics  of  electrons,  ions  and  fluid  to  obtain  spatio-temporal 
profiles  of  electron  density,  ion  density,  electric  potential,  neutral  gas  density  and  neutral  gas  velocity.  We  documented 
how  the  EHD  force  is  controlled  by  the  dielectric  characteristics,  applied  voltage,  frequency  and  the  asymmetric 
configuration  of  the  electrodes,  thickness  of  the  exposed  electrode.  Additionally,  we  utilized  the  force  field  calculated 
from  the  first  principles5  to  predict  the  effect  of  increasing  number  of  electrodes  on  the  wall  jet  while  the  electrodes  are 
kept  at  a  single  phase  and  at  a  phase  gap.  The  momentum  transfer  to  the  gas  was  also  predicted  for  varying  the  distance 
between  the  powered  electrodes. 

Recently  we  have  presented  a  real  gas  model  with  air-like  N2/02  mixture  and  have  shown  preliminary  results 
predicting  the  effect  of  voltage  in  EHD  force  generation  for  eight  species  air  chemistry.11  In  this  paper  we  summarize 
our  recent  developments  with  the  atmospheric  barrier  discharge  model  for  eight  species  of  air-like  medium  and  the 
generated  EHD  force  interactions  over  a  flat  dielectric.  We  also  present  pulsed  dc  and  rf  based  surface  discharge  models 
for  control  of  separated  gas  flow  about  a  conical  cross-section.  This  paper  is  organized  as  follows.  Section  II  describes 
the  problem  statement  for  air  chemistry  and  computed  results.  Section  III  describes  and  interprets  the  flow  control 
simulation  about  the  conical  surface.  Section  IV  summarizes  the  conclusions. 
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II.  Air  Chemistry  in  Plasma  Actuator 

Fig.  1  shows  schematic  of  an  asymmetric  single  dielectric  barrier  plasma  actuator.  It  consists  of  two  electrodes  separated 
by  a  dielectric.  The  upper  electrode  is  exposed  to  the  air  while  the  lower  electrode  is  placed  underneath  the  dielectric. 
The  region  simulated  is  4  cm  long  and  1  cm  high.  The  lower  part  of  the  domain  consists  of  a  0.1  cm  thick  insulator  with 
dielectric  constant  Sd  =  10£0  for  diamond  while  the  upper  part  is  filled  with  air  offy  =  1.0055^0 ,  where  £0  is 

permittivity  of  the  free  space.  We  have  assumed  negligible  thickness  for  electrodes.  The  powered  electrode  extends  from 
x=1.7  cm  to  x=1.9  cm  at  y=0.1  cm,  the  grounded  electrode  is  from  2.1  cm  to  2.3  cm  at  y=0,  with  a  0.2  cm  gap  between 
electrodes  along  the  x-axis.  The  embedded  electrode  is  grounded  and  an  alternating  voltage  of 
(f)  =  1 500  sin(l  0000;zt)  volts  is  applied  to  the  exposed  electrode. 
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m 
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FIG.  1.  Schematic  of  asymmetric  single  dielectric  barrier  plasma  actuators  operating  in  air. 

For  the  air  chemistry  we  neglect  the  metastable  species  along  with  N4+  and  04+  due  to  their  extremely  fast 
recombination  rates.  Also,  the  numerical  complexity  is  further  simplified  by  excluding  nitrous  oxide  at  this  stage.  The 
model  equations  governing  chemistry  of  discharge  are  as  follows.12 


(1)  Nitrogen  model 

(a)  Ionization  e  +  N2  -»  N2  +  2e , 

(b)  Dissociation  e  +  N2  — »  2 N  +  e  , 

(c)  Dissociative  recombination  e  +  N2  —>  N2  , 

(2)  Oxygen  model 

(a)  Ionization  e  +  02  — »  02  +  2e , 

(b)  Dissociation  e  +  02  — »  20  +  e  , 

(c)  Dissociative  attachment  e  +  02  — »  0_  +  O  , 

(d)  Recombination  0~  +  02  ~^0  +  02, 

(e)  Dissociative  recombination  O  +  0~  — >  02  +  e  . 


(la) 

(lb) 

(lc) 

(2a) 

(2b) 
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(2d) 

(2e) 


The  drift-diffusion  form  of  continuity  and  Poisson  equations  for  the  electrons,  ions  and  neutrals  are  as  follows: 
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with  momentum  flux  nava  =  —  sgn(e)na//aV0  —  DaVna ,  and  electric  field  E  =  -V  </> . 

The  electron  temperature  is  calculated  from  E  =  kBTe/(Vne/ne),  which  is  obtained  assuming  an  initial  Boltzmann 
distribution,  ne  oc  exp(e(/> /kBTe).  Various  rate  coefficients  kh  k2,...k%  related  to  in  Eqs.  (la)  through  (2e)  are  obtained 
from  Kossyi  et  al.12  as  functions  of  electron  temperature.  The  secondary  electron  emission  from  the  exposed  surface  is 
taken  as  a  function  of  incident  electron  energy13  and  remains  small,  i.e.,  less  than  10"3.  No  material  sputtering  of  the 
surface  is  considered.  The  current  balance  is  maintained  at  the  surface.  In  Eqns.  (3)-(ll),  na  and  \a  are  density  and 
velocities  of  species  a.  The  mobilities  jua  and  diffusion  rates  Da  are  taken  from  Ellis  et  al.14.  The  bulk  density  of  the  air 
is  taken  to  be  1.3  kg/m3.  Atmospheric  ratio  of  3.6  is  taken  for  nitrogen  to  oxygen  gas  molecules.  The  self-consistent 
formulation  is  solved  using  a  Galerkin  variational  formulation  based  finite-element  method5  to  obtain  electron,  ion  and 
neutral  species  densities  of  nitrogen  and  oxygen,  and  the  electric  potential. 

Initial  and  boundary  conditions  are  as  follows.  The  all  initial  particle  concentrations,  except  those  of  the 
electrons,  nitrogen  and  oxygen  molecules,  are  taken  to  be  zero.  Initial  oxygen  molecule  density  is  taken  as  1024/m3  and 
the  electron  density  is  taken  as  104/m3.  Initial  rf  potential  is  also  zero.  All  neutral  particle  concentrations  are  taken  zero 
at  the  boundaries  of  the  dielectric  and  normal  components  of  flux  are  zero  at  outer  boundaries  of  plasma.  Normal 
components  of  electric  field  are  taken  zero  on  all  the  outer  boundaries  of  the  domain  (dielectric  as  well  as  plasma 
domain).  No  electron  and  ion  current  flow  inside  the  dielectric  is  considered.  Electric  field  normal  to  the  dielectric 
surface  is  discontinuous  by  the  separated  charge. 

The  simulation  results  at  the  quasi-periodic  40th  cycle  are  given  in  Figures  2-11.  Figure  2  plots  the  temporal 
average  of  the  02+  ions  in  the  domain.  The  time  average  was  calculated  by  taking  an  average  of  the  computed  nodewise 
distributed  data  from  21  equally  spaced  time  stations  in  a  2n  cycle.  A  peak  of  2xl013  ions  is  observed  very  close  to 
surface  over  the  grounded  electrode.  The  exposed  electrode  accumulates  an  average  of  1.5xl013  ions  /m3.  The  density 
decreases  sharply  by  two  orders  of  magnitude.  The  thickness  of  this  patch  is  within  a  few  microns  indicating  a  thin 
sheath  layer  (~A,D).  Corresponding  positive  nitrogen  ions  N2+  in  figure  3  show  a  high  concentration  (3xl014  ions/m3) 
very  close  to  the  electrodes.  Interestingly,  the  negative  oxygen  atom  production  shown  in  figure  4  is  of  the  same  level  of 
that  of  positive  nitrogen  molecules.  Obviously,  the  negative  oxygen  atom  plays  a  dominant  role  in  oxygen  chemistry. 
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Figure  2:  Distribution  of  positive  oxygen  ions  shows  a  very  thin  sheath  layer  about  the  electrodes. 
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Figure  4:  Distribution  of  negative  oxygen  atoms  shows  production  of  ions  about  a  very  thin  sheath  (~^D)  layer. 


While  the  nitrogen  and  oxygen  ions  are  concentrated  in  the  sheath  region  about  the  electrodes,  the  neutral 
nitrogen  and  oxygen  atom  production  shown  in  figures  5  and  6,  respectively,  follows  the  electric  field.  The  highest 
density  of  these  neutral  atoms  is  near  the  corner  of  the  exposed  electrode  and  between  the  two  electrodes.  The  oxygen 
atom  number  density  is  four  times  (2x10 14  /m3)  than  that  of  the  nitrogen  atoms.  Figure  7  plots  the  distribution  of  time 


averaged  space  charge  density  («  .  +  n 


-  fl  )  of  the  computational  domain.  Clearly  there  is  a  2mm  thick 


cloud  of  net  positive  charge  of  1013  /m3  above  the  exposed  electrode  with  a  Debye  length  (XD)  thick  negative  layer 
attached  to  the  right  comer.  The  top  of  the  grounded  electrode  is  covered  with  a  cloud  of  negative  charge  of  peak  density 
3.4xlOr 
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/m  .  This  is  reflective  of  a  DC-offset  that  should  be  experimentally  validated  (see  e.g,  Fig.  10  of  Font  et  al.  ). 
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Figure  5:  Time  average  of  nitrogen  atoms  above  show  peak  just  downstream  of  the  exposed  electrode. 
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Figure  6:  Time  average  of  oxygen  atom  production  above  show  peak  just  downstream  of  the  exposed  electrode  is 

nearly  four  times  than  that  of  nitrogen. 
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Figure  7:  Time  average  of  space  charge  density  shows  net  accumulation  of  positive  and  negative  charge. 


The  spatial  variations  of  positive  nitrogen  and  oxygen  ions  along  10  x-locations  are  described  in  Figures  8  (a) 
and  8  (b),  respectively.  The  locations  are  selected  from  1mm  upstream  of  the  exposed  electrode  (x=16mm)  to  2mm 
downstream  of  the  grounded  electrode  (x=25mm).  At  nearly  .5mm  above  the  midpoint  of  the  exposed  electrode,  both 
positive  ions  show  a  peak.  Here  again,  confirming  the  observation  of  Figures  2  and  3,  the  nitrogen  positive  ions  show  an 
order  of  magnitude  higher  number  density  than  the  oxygen  counterpart.  The  sheath  stmcture  is  clearly  visible  in  both 
figures  8  (a)  and  8  (b).  Figure  9  plots  the  time  evolution  of  the  charge  distribution  at  1 1  stations  equally  spaced  in  a  2n 
cycle.  The  selected  x-location  in  this  case  is  at  the  right  comer  of  the  exposed  electrode  (x=19mm).  As  mentioned 
before,  the  oxygen  chemistry  is  largely  dominated  by  the  negative  atoms  that  can  be  clearly  observed  comparing  figures 
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9  (a)  and  9  (b).  The  net  space  charge  density  ( n 


+  n.  —  «  _  —  72  )  due  to  the  species  chemistry  considered  is  shown 

i  t  2  C/2  C/ 

in  figure  9  (c).  While  it  appears  that  the  line  curve  looks  symmetric,  a  close  observation  reveals  a  slightly  negative  bias 
(~1013  /m3)  during  the  cycle.  The  streamwise  (x)  component  of  the  time  averaged  electrodynamic  volume  force 

e(  nN+  +  nQ+  —  nQ_  —  ne  )EX  as  presented  in  figure  10  shows  that  over  the  majority  of  the  actuator,  over  the  exposed  and 

part  of  the  grounded  electrodes,  the  force  is  positive  and  will  push  the  bulk  air  in  the  forward  direction.  Downstream  of 
this,  however,  there  is  a  negative  force  resisting  the  acceleration  of  the  flow.  While  the  accelerating  force  is  about  2100 
N/m3,  the  decelerating  force  is  only  a  fraction  (-600  N/m3).  In  the  experiment  this  may  be  reflected  as  a  dip  in  the  ever 
increasing  velocity  curve  as  time  evolves  though  a  cycle.  We  believe  this  deceleration  present  in  any  dbd  actuator 
inherently  limits  the  performance  and  needs  to  be  carefully  explored. 
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(a)  Population  of  N2+  (b)  Population  of  02+ 

Figure  8:  The  distribution  of  positive  ions  away  from  the  surface  at  2tc/5. 


Figure  9:  Time  evolution  of  positive  and  negative  oxygen  ions  and  the  space  charge  density  at  x=19mm. 


Figure  10:  Time  average  of  streamwise  EHD  force  distribution  based  on  air  chemistry  at  the  40th  cycle. 
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The  distribution  of  computed  electron  temperature  above  the  dielectric  surface  is  plotted  in  figure  1 1 .  It  shows 
the  region  of  high  temperature  at  the  edge  of  the  exposed  electrode  extended  towards  the  grounded  electrode.  This  is 
important  as  the  visible  optical  glow  will  coincide  with  this  temperature  validating  the  experiments.1,2  The  peak  electron 
temperature  is  about  2.2  eV  just  at  the  comer  of  the  powered  electrode  resembling  the  highest  electric  field  in  that 


region. 
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Figure  11:  Distribution  of  electron  temperature  above  the  dielectric  surface  shows  the  peak  temperature  (2.2  eV) 

is  just  downstream  of  the  exposed  electrode. 


III.  Separated  Flow  Control  about  a  Conical  Cross-section 

Schematic  of  flow  control  on  a  fore-body  of  an  aircraft  using  plasma  discharge  plasma  with  an  incident  gas  flow  angle 
a  is  shown  in  figure  12a.  The  geometry  is  obtained  from  Ref.  15.  Width  of  the  domain  is  2  meter  and  height  is  1  meter. 
Gas  is  incident  at  an  angle  a  with  10  m/s  velocity.  Angle  at  the  nose  of  the  forebody  is  5°.  We  have  studied  two 
geometries,  one  with  dielectric  barrier  discharge  plasma  actuator  and  another  with  seven  pairs  of  electrodes  powered 
with  pulsed  dc  voltage.  Thickness  of  electrodes  is  assumed  infinitesimally  small  for  both  the  geometries.  Following  is 
description  of  the  geometry  with  dielectric  barrier  discharge  plasma  actuator.  Dielectric  is  located  from  x=1.2  m  to 
x=1.21  m  with  dielectric  thickness  equal  to  1.0  mm.  The  rf  electrode  is  located  at  x=1.203  m  to  x=1.204  m  at  y=0.301 
m.  Grounded  electrode  is  x=1.206  m  to  x=1.207  m  at  y=0.3  m.  The  electrodes  for  pulsed  dc  geometry  are  located  at 
y=0.3  m  from  x=1.05  to  1.172.  Length  of  each  electrode  is  2  mm  and  distance  between  electrodes  is  8  mm.  There  are 
seven  pairs  of  electrodes.  Each  pair  is  powered  by  a  pulsed  dc  potential. 

Initial  electron  and  ion  densities  are  109/m3.  Initial  voltage  is  zero.  A  10  m/s  neutral  gas  flow  is  imposed  at  the 
left  boundary  with  an  angle  of  attack  of  a  =  15°  (17.5°  with  the  centerline  of  the  forebody).  Initial  air  density  is  1.2 
Kg/m3  and  initial  gas  temperature  is  300  K.  Boundary  conditions  for  Poisson’s  equation  are  as  follows:  potential  <|)  =  4>0 
sin(500007rt)  with  (J)0  =  500  volts  for  DBD  case  and  4>  =  4>0  sin20  (IOOOtcO,  with  (]>o  =  1200  volts  for  pulsed  dc  voltage  case. 
The  potential  is  applied  to  the  first  electrode  of  the  pair  for  pulsed  dc  case.  Other  electrodes  of  each  pair  are  grounded. 
Electric  insulation  condition  (normal  component  of  electric  field  equal  to  zero)  is  assumed  at  outer  boundaries  of  the 
domain  and  at  the  surface  of  forebody.  Following  are  boundary  conditions  related  to  electrons  and  ions  continuity 
equations:  the  currents  flow  normal  only  to  the  electrodes  because  it  is  a  equipotential  surface.  Homogeneous  Neumann 
conditions  are  applied  to  the  outer  edges  of  the  domain  and  electric  insulation  is  assumed  at  the  surface  of  forebody.  The 
currents  flow  normal  as  well  as  parallel  to  the  dielectric  surface.  The  boundary  conditions  for  Navier-Stokes  equations 
are  as  follows:  The  no-slip  condition  (u=0  and  v=0)  is  assumed  for  the  gas  neutrals  at  the  surface  of  forebody,  dielectric 
surface  (for  dielectric  barrier  discharge  case)  and  at  electrodes.  The  velocity  at  the  left  boundary  of  upper  domain  is  set 
equal  to  u=9.65  m/s  and  v=2.58  m/s.  Neutral  boundary  conditions  are  applied  to  the  outer  edges  of  the  domain.  The  self- 
consistent  formulation  is  solved  using  a  Galerkin  variational  formulation  based  finite-element  method  to  obtain  electron 
and  ion  density,  electric  potential,  neutral  velocity,  density  and  gas  temperature. 

Figure  12  (b)  shows  streamlines  and  vectors  of  gas  velocity  at  the  end  of  0.1  sec  without  any  plasma  discharge 
obtained  by  solving  Navier-Stokes  equations.  We  can  see  a  recirculation  bubble  near  the  nose  of  the  forebody.  These 
kinds  of  separation  bubbles  are  harmful  to  the  performance  of  the  aircraft.  An  rf  voltage  for  DBD  or  a  pulsed  dc  voltage 
is  switched  on  after  0.1  second  to  see  the  effect.  Figure  12  (c)  shows  profile  of  pulsed  dc  voltage  as  a  function  of  time. 
The  results  of  Figure  13  are  for  dielectric  barrier  discharge  after  25  cycles  of  plasma  discharge  activity  and  that  of 
Figures  14  for  pulsed  dc  plasma  discharge  after  10  cycles  of  plasma  discharge  activity. 
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(a)  (b)  (c) 

Figure  12:  (a)  Schematic  of  the  conical  section,  (b)  separated  flow  at  17.5  deg  AO  A,  (c)  DC  pulsed  at  1  kHz  rate. 

Figure  13  shows  effect  of  Dielectric  Barrier  Discharge  (DBD)  plasma  actuator  on  flow  control  on  a  forebody  of  an 
aircraft.  We  do  not  solve  temperature  equations  (3)  for  the  results  of  this  figure.  After  0.1  second  rf  voltage  §  =  500 
sin(27i ft)  volts  and  /  =  25  kHz,  is  switched  on  to  see  the  effect.  Figure  13  (a)  shows  vectors  of  time  averaged 
electrodynamic  force  e(wf  -  ne) E.  It  can  be  seen  from  the  figure  that  force  is  directed  in  forward  direction  which  may 
help  in  elimination  of  recirculation.  Figure  13  (b)  shows  streamlines  and  vectors  of  gas  velocity  after  25  cycles  of  rf 
discharge  activity.  We  can  see  that  recirculation  is  affected  by  DBD  plasma  actuator;  however,  flow  is  not  attached  to 
the  surface  all  over  the  surface  of  forebody.  This  is  because  single  dielectric  barrier  plasma  actuator  is  not  sufficient  to 
cover  the  entire  length  of  recirculation  and  recirculation  is  not  eliminated.  Multiple  dielectric  barrier  discharge  plasma 
actuators  on  the  surface  of  forebody  may  control  vortex  dynamics  of  the  flow.  We  suggest  multiple  pair  of  electrodes 
powered  with  pulsed  dc  potential  as  an  alternative  to  multiple  DBD  plasma  actuators  in  this  paper. 


(a)  (b) 

Figure  13:  (a)  Time  averaged  EHD  force  vectors  shows  force  distribution  about  the  DBD  actuator,  and  (b)  Gas 
flow  streamtraces  after  50  cycles  of  barrier  discharge,  (x  and  y  are  in  meters) 

Figure  14  shows  the  effect  of  plasma  discharge  on  flow  control  about  a  the  crosssection  of  a  conical  forebody 
due  to  electrodynamic  force  only.  We  want  to  see  the  effect  played  by  electrodynamic  force  on  flow  control  on  a 
forebody.  Figure  14  (a)  shows  charge  separation  and  vectors  of  electrodynamic  force  at  the  peak  of  the  spike.  At  its  peak 
the  EHD  force  vectors  are  directed  mostly  downwards  which  helps  in  destruction  of  the  big  recirculation  bubble  and 
formation  of  small  bubbles.  Figure  14  (b)  shows  streamlines  and  vectors  of  gas  velocity  after  10  cycles  of  plasma 
discharge  activity  by  a  pulsed  dc  potential.  We  can  see  that  effect  of  plasma  discharge  on  recirculation  is  highest  in  this 
case  and  a  big  recirculation  bubble  is  destroyed  and  small  bubbles  are  formed.  These  small  bubbles  may  not  be  that 
much  harmful  to  the  performance  of  the  aircraft  as  a  big  bubble. 
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Figure  14:  (a)  Space  charge  distribution  at  the  peak  of  pulsed  dc,  (b)  Streamtraces  show  the  effect  of  EHD  force 
on  flow  structure. 


IV.  Conclusions 

We  have  utilized  a  finite  element  based  formulation  to  model  the  air  chemistry  around  a  plasma  actuator  and 
numerically  compared  the  capability  of  a  set  of  pulsed  dc  and  a  single  dbd  actuator  to  control  largely  separated  flow 
about  a  forebody  conical  cross-section.  The  N2/02  air  chemistry  consists  of  eight  charged  and  neutral  species  and 
demonstrated  the  dominance  of  N2+  and  O'  ions.  The  time  average  of  the  streamwise  force  shows  mostly  acceleration 
above  the  actuator  but  a  small  decelerating  force  downstream  of  the  powered  electrode.  This  will  induce  a  fluctuation  in 
the  temporal  evolution  of  streamwise  velocity.  This  decelerating  component  needs  further  investigation  as  this  may 
inherently  limit  the  performance  of  a  dbd  plasma  actuator.  The  computed  electron  temperature  based  on  the  electric  field 
resembles  the  visible  glow  that  is  commonly  reported  in  the  literature.  For  the  conical  cross-section  of  a  forebody,  a 
large  separation  occurs  due  to  a  10  m/s  at  17.5  deg  AOA  neutral  gas  flow  past  the  geometry.  A  single  dbd  actuator  could 
not  control  such  separation.  Application  of  seven  pulsed  dc  actuators  worked  significantly  better  to  control  the  large 
vertical  structures.  In  near  future  we  will  investigate  whether  this  effect  is  due  to  dc  Joule  heating  of  the  boundary  layer 
(cf.  Maslov  et  al.15),  electric  force  or  a  combination  of  them. 
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